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Abstract. The OGLE database is revisited to investigate in more detail the properties of the Fourier parameters. 
Methodological improvements led us to identify a clear separation among High-Amplitude S Scuti (HADS), RRc 
and RRab stars. The bimodal distribution of the _R2i parameter in If ADS stars is explained as a contamination 
effect from RRc stars: there is evidence that all stars with 0.20< P <0.25 d are RRc variables. The previously 
claimed existence of a subclass of unusual HADS is demonstrated to be a spurious result. Candidate overtone 
pulsators are found among HADS and RRc variables. The properties of the Fourier parameters are discussed as a 
function of the physical conditions in the stars involved. Among the field RRab stars we detected different light- 
curve groups producing distinct "tails" in the Fourier plots for P >0.55 d; evolutionary phases or the combination 
of different physical conditions (not only metallicity) are suggested to explain this separation, observed also in 
the cluster RRab stars. The stellar parameters of RRc stars in a given globular cluster show different tendencies 
than those of RRc stars from different clusters. 

Key words. Methods: data analysis - Techniques: photometric - Astronomical data bases: miscellaneous - Stars: 
oscillations - Stars: variables: S Set - Stars: variables: RR Lyr 



1. Introduction The Fourier decomposition is a powerful method to 

" " ' describe in a synthetic and safe way the shape of a light 

The Optical Gravitational Lensmg Experiment (OGLE; consists in fitting the photometric measurements 

Udalski et al. 1993) is a project whose main goal is a j^jgg^jj^g Qf series 
systematic search for microlensing events. As a valuable 

by-product, a large amount of photometric data was ac- n 

quired on a variety of variable stars. As demonstrated by I{t) — Aq + ^ Ai cos[27ri(t — To)f + (f)i] (1) 

several studies, the light curves of pulsating stars depend i=i 

on their physical parameters. Moreover, the observed light 

1 J .,1 ,1 J • J f u J J where I(t) is the magnitude observed at time t, An the 

curves can be compared with those derived from hydrody- ^ ' , ° . ' " 

1 ; • , J , r u 1.1 mean magnitude, Aj the amplitude of the i-component 

namic pulsational models, iheretore, once an observable t> i « f f 

quantity has been recognized as depending on stellar pa- ' " ^ harmonic), / the frequency (/=1/P, where P 

rameters, it is possible and useful to fix the relation be- '^J^"" P^"^^ ^^S^* variation), the phase of the 

tween this quantity and the light curve. That can help us ^-component at t - To. 

in obtaining stellar parameters for stars for which we have ^he Fourier parameters can be subdivided into two 

only hght curves, i.e. the by-product of the microlensing groups: the amplitude ratios i?y = A,/A^ (i.e. i?2i = 

projects. M/Ai, R31 = A3/A1, R32 = A3/A2, ...) and the phase dif- 
ferences 4)j.i = i(t)j -j(t>i (i.e. (t)2i=(t>2 - 2</>i, 031 =(^3 - 301, 

Send offprint requests to: E. Poretti '^32 = 203 - 302, ..•)■ Note that the use of the sin term 

* Tables 3, 4 and 5 are only available in elec- mstead of the cos leads to results that are not directly 

tronic form at the CDS via anonymous ftp comparable, owing to the it/2 shift of each component. 

to cdsarc.u-strasbg.fr (130.79.128.5) or y ia Because it provides quantitative parameters to define the 
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a suitable tool for classification purposes, once care is 
taken to avoid spurious results or misinterpretations of 
features of the light curves. Moreover, the plots of the 
Fourier parameters versus periods showed the presence of 
discontinuities in the progressions, which can be ascribed 
to resonances between different pulsational modes, as in 
the case of Classical Cepheids at P~10 d. These diagrams 
are therefore useful tools to perform asteroseismology in 
giant stars, as the Cepheids (Poretti 2000a), even if some 
caution is necessary to locate exactly the position of the 
resonances (see Kienzle et al. 1999 and Feuchtinger et al. 
2000 for a discussion of the resonance between the fourth 
and first overtones in the case of s~Cepheids). 

In this paper, we examine the Fourier plots of the pul- 
sating stars with period less than 1 d contained in the 
OGLE database in order to emphasize their properties. 

2. The OGLE database 

We considered the 268 pulsating stars included in the 15 
fields covered by the OGLE database. The number of stars 
is large enough to constitute a conspicuous sample, able to 
provide unambiguous characterizations of these variables, 
and small enough to allow a star- by-star analysis, thus 
avoiding the misinterpretations that could arise from au- 
tomated analysis. Each light curve is suitable to be Fourier 
decomposed thanks to the high number of homogeneous 
measurements (with a few exceptions) and we expected 
a clear detection of features and peculiarities, both for 
individual stars and for classes of variables. The analy- 
sis consisted of the separate steps described in the next 
subsections. 

2.1. The verification of previous results 

The primitive idea was simply to investigate the Fourier 
plots obtained by Morgan et al. (1998a, 1998b, 1998c) to 
clarify some ambiguous patterns. In particular, the 031— P 
plot does not show any clear progression; the (j)2i—P plot 
also shows a strange scatter at short periods. If these re- 
sults were confirmed, they should be considered as a fail- 
ure of the capabilities offered by the Fourier decomposition 
as a sharp tool to describe light curves. 

A close examination of the coefficients of the least- 
squares fit reported by Morgan et al. (1998a, 1998b, 1998c) 
reveals that in many cases the small amplitude terms 
are not significant, the error bars being as large as the 
amplitudes themselves. In other cases, the presence of 
large-amplitude high-order harmonics (i.e. Ai^i/Ai ~ 1) 
implies a really odd fitting light curve with humps and 
bumps. In general, Ai^i/Ai < 0.5 should be the rule far 
from resonances. Only the necessity of fitting the double 
maximum of RRc variables, the well-defined humps at 
minimum light of some RR Lyr variables or very asym- 
metrical curves should involve deviations from this rule, 
generating a few consecutive terms having similar ampli- 
tudes. Therefore, we matured the idea that a re-analysis 
of all the datasets was necessary to handle more reliable 



Fourier coefficients. At this point, we were forced to scru- 
tinize each star again. 

2.2. The fit of the OGLE light curves 

We excluded 6 stars from the original sample (BW1_V31, 
BW9.V55, BWC.V22, BWC.V28, BWC.V97 and 
BWC_V56). For these stars, the number of measurements 
is less than 60 and the fit is quite uncertain, hampering 
a reliable determination of the Fourier parameters. In 8 
other cases the measurements do not show significant de- 
viations from a perfect sine-shaped light curve, probably 
owing to the small amplitude: BW1_V90 (full-amplitude 
0.08 mag), BW2.V44 (0.10 mag), BW4.V55 (0.04 mag), 
BW4.V94 (0.06 mag), BW5.V117 (0.10 mag), BW6_V120 
(0.10 mag), BW7_V74 (0.05 mag) and MM5-B.V141. It 
should be noted that 0.049 d< P <0.180 d for aU these 
stars (and for 6 of them P < 0.10 d): the shorter periods 
are compatible with the possibility that these stars are 
overtone pulsators. Figure |^ shows a light curve in which 
it was not possible to detect the amplitude and phase 
values of the first harmonic in a reliable way. Morgan et 
al. (1998c) considered four of these variables (BW1_V90, 
BW4.V94, BW6.V120, BW7.V74; see their Tab. 1) as 
unusual S Set. This misleading conclusion was due to 
the over-evaluation of the significance of their 02i, 03i, 
041 parameters. 

In 8 more cases (BW1.V7, BW4.V132, BW5.V135, 
BW5.V174, BW9.V35, BW10.V39, BWC.V65 and MM7- 
B_V7) it was not possible to apply the Fourier decompo- 
sition: a high-order fit was requested to interpolate the 
measurements describing a very asymmetric shape, but 
the result is an improbable curve with evident bumps. 
In such cases the combination of the noise with the large 
number of free parameters is responsible for the unaccept- 
able (from an astrophysical point of view) curve, even if 
the least-squares routine brings out figures. The frequency 
analysis did not reveal any second periodicity. Figure 
shows an example. 

Finally, we reduced to 246 the number of stars whose 
light curves were able to be decomposed. The 22 stars we 
rejected can be considered as partially responsible for the 
unclear Fourier plots obtained by Morgan et al. (1998a, 
1998b, 1998c). In particular, this led us to a first spe- 
cific result. We had already found that 4 of the stars de- 
scribed by Morgan et al. (1998c) as unusual d Set stars 
have a sine-shaped light curve. Proceeding further, we 
found that the curves of all the remaining "unusual" stars 
could be fitted by a M=2 Fourier series and hence the 
031 and 041 values don't have any significance. Since the 
02iparameters of these remaining stars assume quite nor- 
mal values, the claimed subclass of unusual 6 Set stars 
previously reported doesn't exist. 

Moreover, we changed the decomposition order M of 
about 70% of stars; in this way the Fourier plots involving 
high-order terms improved in clearness. 
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Fig. 1. Example of a sine shaped light curve (BW4_V55, 
P=0. 15958 d). The small amplitude (the fitting curve has 
a full amplitude of 0.04 mag) and the noise (rms residual 
0.016 mag) do not allow the detection of a significant first 
harmonic contribution. 



17.4 
I 

17.6 
17.8 




0.2 










0.5 
Phase 


1 




1 ' 


1 1 1 1 1 


1 1 1 1 1 1 1 
BW5_V135 


1 1 1 


' 1 




















1 . . . , 1 . . . . 1 , , . , 1 . . , . 1 





5 


10 

Freq. [c/d] 


15 


20 



Fig. 2. Top panel: example of an unacceptable least- 
squares fit (BW5_V135, P=0.58825 d). It is hard to be- 
lieve that the bumps are real; the asymmetry of the light 
curve requires a high order of the Fourier fit (M =4) and 
the large scatter enhances the amplitude of the high-order 
terms. Bottom panel: the power spectrum of the residuals 
does not show any significant second periodicity. 



2.3. The period refinement 

The Fourier decomposition of the light curves by a sin- 
gle frequency (plus harmonics) could not always give a 
satisfactory fit, sometimes leaving a high rms residual. In 
such cases, it is quite obvious to investigate the datasets 
searching for other periodicities; this examination had not 
been done in the previous analysis. In many cases our fre- 
quency analysis revealed the presence of other peaks in 
the power spectrum. Since most of these peaks arc ob- 
served very close to the main one and since such periods 
are expected as signatures both of nonradial oscillations 
and of a Blazhko effect, we must be sure that the reported 
periods are accurate. If not, a wrong period can generate 
spurious alias peaks or change the observed value of the 
neighbouring real peaks. 

Therefore, we decided to refine the periods of all the 
stars and then to re-analyzc in frequency all the datasets. 
To this end, we considered as preliminary solutions the 
least-squares fit obtained by using the period values re- 
ported in the OGLE database and the order M established 
in the previous step. This solution was given as input pa- 
rameters to a code keeping locked the relations between 
the main frequency and the harmonic terms (MTRAP; 
Carpino et al. 1987); the best fit was searched for around 
the preliminary solution. Once the correct period was 
obtained, we re-examined in frequency all the datasets 
to detect secondary periodicities. No prewhitening was 
done: only the frequency values of the main term and its 
harmonics were considered as input values (known con- 
stituents), not their amplitude and phase. This approach 
is the same as the one used by Pardo & Poretti (1997) to 
detect the frequency content of double mode Cepheids. 
The period refinement was marginal in most cases: for 73 
stars (30% of the sample) it is smaller than 1 • 10~^ d, for 
168 stars (72% of the sample) is smaller than 5 • 10~^ d. 
The mean error on the period is 1.3 • 10~^ d and for about 
100 stars out of 234 the refinement was made within this 
error bar. As a by-product, wc made a further revision of 
the significance of the small amplitude terms. The main 
result of this search is described in the next subsection. 



2.4. Ttie different frequency content of pulsating 

variables 

At the end of this step we separated the monoperiodic 
variables from the other stars. Of course, a second, small- 
amplitude term may always be hidden in the noise, even 
if the fit with a single-period looks satisfactory. However, 
after the further check we made, we can be confident that 
its effect on the main period should be negligible. As a 
result of the period refinement procedure and the subse- 
quent frequency analysis, we evidenced 53 stars showing 
a well-established multimode behaviour: 

1. double-mode stars pulsating in the fundamental radial 

mode and in another radial overtone (7 candidates), 
with a well-separated pair of frequencies; 
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Table 1. Averaged errors on the 02i, ' 
rameters. 
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Fig. 3. Histogram of the standard deviations of the reh- 
able Fourier fits of the hght curves of the 192 monoperiodic 
stars. 



2. stars showing one single peak close to the main fre- 
quency (30 candidates), with a ratio /i 7/2^^0. 99 or 
~1.01; 

3. stars showing two peaks symmetric with respect to the 
main frequency and, in some cases, to the harmonics 
(7 candidates); 

4. stars showing a second peak very close to the main 
one (9 candidates), /i//2'^0.999 or '^l.OOl, suggesting 
a period variation or an unresolved, very long Blazhko 
effect. 

All these types of multimode variables will be discussed in 
forthcoming papers. Not considering the Pop. II Cepheid 
BWC_V1 (P= 1.74797 d), we got a rehable Fourier de- 
composition of 192 monoperiodic variables with P <1 d. 
In general, the fit is satisfactory even if some light curves 
appeared to be noisy. That can be ascribed both to dif- 
ficult image reduction (close companion, nearby bright 
stars, ...) or to slow instrumental drifts, which appear in 
the frequency analysis as peaks very close to 0.000 cd~-^ or 
1.000 cd~^. In 136 cases (i.e. the 71%) the least-squares 
fits have a standard deviation better than 0.03 mag 
(Fig. a value quite acceptable for stars of /-magnitude 
15-17 measured on a time baseline of three years. 

We would like to emphasize the fact that all the work 
described above constitutes a complete independent re- 
analysis of the whole OGLE database. The new multi- 
periodic pulsators we found, the sharper revision of the 
light curve fits and the elimination of unclear subclasses 
of variables configure our Fourier analysis of the OGLE 
database as a very reliable tool to investigate the proper- 
ties of the pulsating variables with P < 1 d. 
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3. The Fourier diagrams: the separations between 
variable classes 

We can finally handle the well-defined Fourier plots shown 
in Figs. ^, H, H and |^. Table [l] lists the averaged errors on 
the Fourier parameters for each class, as calculated from 
the errors on the least -squares fits. 

Considering P <1 d, the pulsating stars should be sub- 
divided into High- Amplitude S Scuti (HADS), RRc and 
RRab stars. However, the continuity of the progressions of 
the phase parameters (Figs. ^, ||, ^ masks the separation 
between HADS and RRc stars. Such a subdivision is more 
appreciable in the R21—P plot (Fig.0). As a first step, a 
close examination can help in separating variable classes 
more precisely. For the sake of clarity, the subdivisions 
are shown right away in Figs. Ill and I by means of 
different symbols. 

3.1. HADS or RRc ? 

In the period range 0.15-0.25 d, separating an RRc vari- 
able from a HADS is not an easy task, since for example, 
the 021— P plot (Fig. I shows a merging of two fami- 
lies of points, without any apparent jump. However, the 
light curve of a star following the progression coming from 
the shortest period shows an asymmetrical light curve 
with a sharp maximum (BW6_V30, P=0. 18902 d; Fig. | 
top panel), while in the opposite direction we generally 
find more smoothed light curves with rounded maxima 
(BW9.V114, P=0.20767 d. Fig. | bottom panel). If the 
difference is not quite evident in the 02i— P and 03i— P 
plots (Figs. ^ and ||), it can be better appreciated both in 
the 041 -P plot (Fig. |) and in the R21-P one (Fig. 0). 

We didn't detect any A4 term for stars in the 0.20 d< 
P <0.25 d interval, suggesting that the light curves are 
more sinusoidal: it looks natural to ascribe this feature to 
the rapid evolution of the light curves of RRc stars (see 
the steady slope of 04iparameters for P > 0.30 d). In the 
latter, a well defined dip is observed in the interval from 
0.20 to 0.40 d, suggesting that all these stars belong to 
the same family. In particular, note the P2i>0.2 values 
for P <0.20 d and the P2i<0.15 for 0.20 d< P <0.25 d. 
We can conclude that the Fourier parameters of stars in 
the interval 0.20 d< P <0.25 d are naturally linked with 
the longer periods (i.e. RRc stars) rather than with the 
shorter ones (i.e. HADS stars). 
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Fig. 4. Fourier parameter (/)2iagainst Period for the 192 
monoperiodic stars with P <1 d in the OGLE database. 
Filled circles: HADS; open triangles: RRc; crosses: RRab 
stars. 
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Fig. 6. Fourier parameter (/)4iagainst Period for the 192 
monoperiodic stars with P <1 d in the OGLE database. 
Same symbols as in Fig. ^ 
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Fig. 5. Fourier parameter (/)3iagainst Period for the 192 
monoperiodic stars with P <1 d in the OGLE database. 
Same symbols as in Fig. ^. 
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Fig. 7. Fourier parameter P2iagahist Period for the 192 
monoperiodic stars with P <1 d in the OGLE database. 
This plot yields the clearest separation between HADS 
and RRc stars. Same symbols as in Fig. 0. 
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Fig. 8. Top panel: note the sharp maximum in the 
hght curve of BW6_V30 (P=0. 18904 d), which can be 
considered the HADS with the longest period. Bottom 
panel: note the rounded maximum in the light curve of 
BW9_V114 (P=0.20767 d), which seems to be the extreme 
short-period RRc. 



Another way to separate HADS from RRc stars is 
to plot the Af term versus the A2 one (Fig. ||, see also 
Antonello 2000). There is some crowding close to the ori- 
gin, but the two classes are well separated in the critical 
range 0.15-0.25 d. 

Therefore, we have demonstrated that when dealing 
with a well-defined sample of accurate light curves, it is 
possible to separate HADS from RRc stars not on the basis 
of a mere assumption (Antonello 2000), but by looking 
at the amplitude values. We shall find below a further 
confirmation of this separation. 

The separation between HADS and RRc stars by using 
P=0.20 d as a clean cut cannot be considered as an invi- 
olable rule; as an example, the symmetrical light curve of 
MM5-B_V141 (P=0.181 d) suggests that this star belongs 
to RRc variables rather than to HADS. We also note that 
the presence of second overtone pulsators among RRc vari- 
ables is always an open possibility; such pulsators should 
show more symmetrical light curves and the low R21 val- 
ues observed here can be a hint. 




0.01 0.02 
A,' [mag^] 

Fig. 9. The amplitudes of the first Fourier terms separate 
HADS (circles) from RRc (triangles) stars. Filled circles 
are HADS stars with 0.15< P < 0.20 d; filled triangles 
are RRc stars with 0.20< P < 0.25 d. 

3.2. RRc or RRab ? 

Looking at Fig. ^, RRc and RRab stars are clearly sep- 
arated. In general, RRc stars have a period P <0.35 d 
and for this period (/)2i>4.5 rad; RRab stars have a pe- 
riod P >0.42 d and for this period 02i<4.O rad. Similar 
separations can be found in all the phase and amplitude 
plots. 

There are only a few exceptions, i.e. the presence 
of an isolated short-period RRab star, MM5-B_V31 
(P=0. 36456 d) and of three long-period RRc stars, 
BWl.Vll (P=0.38434 d), BWC.V81 (P=0.38588 d) and 
BW2_V8 (P=0.39402 d). However, the differences in the 
light curves are quite evident also in these intermediate 
cases (Fig. |l^). The R21 — P plot provides the best sepa- 
ration between short -period RRab stars and long-period 
RRc stars (Fig. While the RRab star displays normal 
values for its class, the three RRc stars display intermedi- 
ate values, but always closer to RRc- than to the RRab- 
type (especially in the P2i~P plot). It should be noted 
that there are few stars in the period interval from 0.35 d 
to 0.40 d; therefore it is probable that these stars have 
some particularity. 



4. The Fourier diagrams: 
different classes 



insights into the 



The large number of stars forming each class allows us to 
discuss some aspects in detail. 

4.1. HADS stars 

The 02ivalues are quite constant from 0.07 d to 0.20 d, 
but they show a slight trend to decrease toward shorter 
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Fig. 10. Top panel: MM5-B_V31 is the RRab star hav- 
ing the shortest period (P=0. 36456 d). Bottom panel: 
BW1_V11 is the RRc star having the longer period 
(P=0. 38434 d). Note also the strong difference in ampli- 
tude. 



periods. For the fl stars with P <0.07 d, we get 
02i=3.92±O.O6 rad, while for the 28 stars with 0.07 d< 
P <0.20 d we get 02i=4.23±O.O4. This tendency was al- 
ready noticed in the uj Cen stars (Poretti 1999), where 
the presence of stars with shorter period (down to 0.036 d) 
makes the trend more evident. When considering the vari- 
ables with P <0.05 d, the slope of the (j)2i—P progres- 
sion is not only confirmed, but seems to become steeper 
again, suggesting a second break (see also Fig. 1 in Poretti 
2000b). The change in the slope can also be seen in the 
4>zi—P and (j)4i—P plots, but it is difficult to place it. The 
gap in the period between 0.080 d and 0.094 d is a further 
complication. It should be noted that such a gap is not 
filled by considering variables we excluded for some reason 
or other: the only star in this interval is the double-mode 
pulsator BW1.V207. 

In correspondence with the break in the 02i~-P pro- 
gression, Antonello (2000) noted a saddle (centered at 
P=0.10 d and flanked by two maxima at 0.08 d and 0.14 d) 
in the amplitude plots. He correlated this feature with the 
stage of stellar contraction before the expansion of the 
envelope. As a matter of fact, the change of slope is un- 
likely to be due to a resonance effect or a different pulsa- 
tion mode, which both should produce much more evident 




Fig. 11. Histogram of the i?2ivalues observed in the light 
curves of HADS stars. The distribution does not show any 
peculiarity considering stars with 0.20 d< P <0.30 d as 
RRc stars. 



changes in the light curves. As regards the latter possi- 
bility, we remind the reader of the presence of perfectly 
sine-shaped light curves in this period range, which could 
be the signature of an overtone radial mode. 

Another feature of interest is the bimodal distribution 
of the i?2i values found by Antonello et al. (1986), ques- 
tioned by Hintz & Joner (1997), re-confirmed by Musazzi 
et al. (1999), simultaneously weakened by Morgan et al. 
(1998c) and strengthened by Templeton et al. (1998), 
to be finally found unclear by Templeton (2000). When 
plotting the values obtained from the OGLE stars with 
P <0.30 d (as in the Antonello et al. paper), we observe 
a bimodal distribution (with different maxima). However, 
when considering the stars with P <0.20 d (i.e. the gen- 
uine HADS stars, as previously established), the distri- 
bution loses that characteristic, showing a single maxi- 
mum toward high R21 values (Fig. [Tl| ). Such high val- 
ues are required to explain the strongly asymmetric light 
curves of HADS stars. Therefore, the bimodal distribu- 
tion is generated by the low i?2 lvalues typical for RRc 
stars having short periods, as seen above. This is also con- 
firmed by looking at the Tab. 2 reported by Antonello et 
al. (1986): the i?2i<0.20 values are mostly found in stars 
with P >0.20 d. As an interesting feedback, this results 
strengthens our confidence in considering i?2i=0.20 as the 
separation between HADS and RRc for periods around 
0.20 d. 
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4.2. RRc stars 

There are no peculiarities in the Ught curves of RRc stars. 
The only exceptions are the short- and long-period stars 
which are respectively close to HADS and RRab variables 
and whose parameters show intermediate values in some 
cases. The Fourier parameters are well confined in narrow 
loci in all the plots (Figs. ^ |[ || and |^) and their progres- 
sions describe very well the main characteristic of these 
stars, i.e the slow but continuous change in the shape of 
the light maximum, which appears to be double or flat in 
most cases. When the light curves are more rounded, i.e. 
around 0.31-0.32 d, there is a corresponding dip in the 
^31 — P and Rii — P diagrams. 

4.3. RRab stars 

The RRab class is richer in particularities than the RRc 
one. In the phase diagrams there is an unique, narrow 
progression from 0.40 d to 0.55 d. After this uniform be- 
haviour, we see the (f)ii values diverge into three "tails" . 
The Fourier 04iparameters provide the clearest separation 
among tails (Fig. |l^). The separation is usually larger than 

0. 50 rad, while the fits of the light curves give a mean er- 
ror on the 04iparameter of 0.16, 0.14 and 0.08 rad for the 
upper, central and lower tails, respectively. In a few cases, 
looking at a particular ipji — P plot, it is not easy to de- 
cide which tail a star belongs to, owing to the error on the 
phase parameter we are examinating. However, this task 
can be done in a much more reliable way by considering 
all the 4>ji — P plots; in that way we take full advantage 
of the presence of the three tails in each plot. As proven 
by the three panels in Fig. |l^, the same stars form the 
same type of tail, though with a different degree of clar- 
ity. The upper and central tails are merged together for 
P <0.55 d and cannot be separated; the lower tail origi- 
nates at P ~0.50 d. Moreover, the OGLE sample does not 
contain many stars with P > 0.65 d; therefore, this part 
is not well defined in each tail (see Sects. 5.1 and 5.2). 

Once evidenced in the Fourier plots, the differences 
can be evaluated also by a direct examination of the light 
curves (Fig. |l3|) : 

1. The lower tail is formed by stars having the steepest 
rising branches and well defined bumps before them; 
the max;ima are sharp and the descending branches are 
concave (Fig. ^ lower panels). Moreover, there are 
no great changes when the period increases (compare 
the light curve of BW11_V3, F=0.596 d, with that of 
BW7_V48, P=0.639 d). Stars belonging to this tail can 
be recovered down to P=0.523 d. 

2. The light curves of the stars forming the central 
tail show a less steep ascending branch; the descend- 
ing branch has a sort of corrugated shape (Fig. p3| , 
middle panels; note that the points, not the fitting 
curve, are describing it) for shorter periods (BW9_V14, 
P=0.607 d) and shows a tendency to become linear for 
longer periods (BW4_V12, P=0.639 d). Moreover, the 
bump is wider and less marked than in the lower tail. 



3. The upper tail is formed by light curves having the 
slowest rise to maximum (Fig. upper panels). The 
bump at minimum light is visible at short periods 
(BW1_V25, P=0.600 d), but disappears at longer ones 
(BW8_V20, P=0.677 d); moreover, the descending 
branch looks convex. The amplitude is smaller than 
in the two previous cases and a 2"'^-order fit is often 
sufficient. Consequently, the R21 — P plot clearly sepa- 
rates the stars belonging to this tail from other RRab 
variables; in Fig. |^ they originate the points below 0.4 
after P=0.60 d. 

As the general shape of the curves changes from one 
tail to another, all the phase parameters are sensitive to 
differences, particularly the higher-order ones, since they 
should fit the bump at minimum and the differences in 
skewness. This explains why the separation of the tails is 
more conspicuous in the 041— P and ^51 — P plots than in 
the 021— P and 031— P plots. 



5. Discussion 

The complete re-analysis of the OGLE database gave us 
the opportunity to demonstrate the non-existence of the 
unusual 5 Scuti stars reported by Morgan et al. (1998c), to 
recognize the bimodal distribution of the P21 parameters 
in the HADS sample as a contamination effect, to separate 
in a better way the three classes of variable stars. Other 
aspects deserve further discussion, also considering both 
the inputs needed by theoretical work and the necessity 
of obtaining clear indications from monoperiodic stars in 
order to better understand multimode pulsators. 

5.1. Features in the progressions 

The Fourier parameters are confirmed to be a powerful dis- 
criminator between the first overtone mode of RRc stars 
and the fundamental mode of RRab stars. The same sepa- 
ration is discernible, though with more difficulty, between 
the RRc stars and the HADS stars, considered to be F 
mode pulsators. These latter stars show a changing slope 
at P '^O.IO d; the possibility of a presence of lO pul- 
sators among short-period HADS is suggested by some 
sine-shaped light curves and by a small group of stars 
having high values of 02i (Poretti 1999). 

None of these progressions shows a clear, abrupt 
change suggesting the presence of a resonance. We can 
speculate that the upper tail of the RRab changes at 
P =0.64 d (Fig.|l2|), as there are two stars (MM5-B_V20 
and BW10_V95) showing a P2ivalue and a shape typical 
of upper sequence stars, but a significantly lower 04i value. 
Moreover, the 05i and 06i parameters show a regular pro- 
gression, confirming that the two stars belong to the upper 
tail. It will be interesting to verify the hint of a resonance 
involving the 3rd overtone by means of a larger database. 
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Fig. 12. Fourier ip^i, 04iand 05i parameters of RRab stars 
in the OGLE database (/ light). After a common progres- 
sion mitil P=0.50-0.55 d (filled squares), the 4>ji values 
diverge into different "tails" . The symbols used to identify 
the three tails are the same as those used in Fig. to plot 
the different light curves. 



5.2. The three tails of the RRab progression 

Kovacs & Jurcsik (f996, 1997) and Jurcsik (1998) have 
derived formulae that estimate the absolute magnitudes, 
metallicities, intrinsic colours and temperatures of RRab 
stars from the periods, amplitudes and phases. In partic- 
ular, the [Fe/H] value is related to P and (j)3i. Therefore, 
the differences in the progression naturally transpose into 
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Fig. 13. Light curves of different RRab stars belonging to 
the three different tails visible in the phase Fourier plots. 
The upper panels (squares) show a short- and a long- 
period RRab characterized by high values of the phase 
differences (upper tail) . The central panels (circles) show a 
short- and a long-period RRab characterized by medium 
values of the phase differences (central tail). The lower 
panels (crosses) show a short- and a long-period RRab 
characterized by low values of the phase differences (lower 
tail). The RR Lyr variables in each column have approx- 
imately the same period and thus differences in the light 
curve shape at the same period can also be assessed. Note 
also the differences in the amplitudes. 

a difference in metallicity. After transforming the phase 
differences values from / to light using the equations 
provided by Morgan et al. (1998b), we indeed obtained 
a slightly different mean metallicity: the 20 stars of the 
upper sequence yield [Fe/H]=-0.96±0.04, the 17 stars of 
the central one -1.39±0.05 and the 9 stars of the lower 
one -1.50±0.08. We note that they are strictly monope- 
riodic pulsators, and the condition of regular light curves 
(Kovacs & Kanbur 1998) is then largely satisfied. 

To investigate in more detail the possible effect of the 
metallicity, we revisited the (pn—P progression for the RR 
Lyr variables in NGC 6171 (Clement & Shelton 1997; we 
omitted two unclear cases), NGC 6362 (Olech et al. 2001), 
M3 and M9 (Clement & Shelton 1999), M55 (Olech et al. 
1999), M5 (Kaluzny et al. 2000a) and NGC 6934 (Kaluzny 
et al. 2000b). We transformed the phase values obtained 
from a Fourier series in sin into phase values from a cos 
series. The parameters obtained in this way are plotted in 
Fig.jlJ using the same symbols as in Fig. ^to indicate the 
different tails. Note the sharpness of the OGLE progres- 
sions compared with the progressions relative to cluster 
stars; observations carried out with different instruments 
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Fig. 14. Fourier parameters of RRab stars in some 
globular clusters {V light). The central tail (filled and open 
circles) is composed by stars from NGC 6362, M55, M9, 
M5, NGC 6934 and M3. The upper tail (open squares) is 
composed by all the stars from NGC 6171 and few stars 
from NGC 6362, M5 and M3. The lower tail (crosses) in- 
dicates stars from M55 and NGC 6934. 



and much more crowded fields can explain the larger scat- 
ter. 

Most of these globular clusters show overlapping cen- 
tral tails (NGC 6362, M9, M5 and NGC 6934), despite a 
large spread in metallicities. However, looking at Table H 
we can see that the expected general distribution is re- 
spected: metal-rich globular clusters have some or most 
stars on the upper tail and none on the lower one, while 
the RR Lyr stars of the metal-poor clusters populate the 
lower tail. Note also that the two extreme cases (NGC 
6171, metal-rich, and M55, metal-poor) describe the two 
progressions well above and below the central one, respec- 
tively. 

However, all the stars of the same cluster should have 
the same metallicity and therefore the 04i— P plot of each 
cluster should show a single tail, but in three cases we 
actually observe the co-existence of two kinds of tails (see 
also Tab. 1): 

1. in M3 there is a tendency for long-period RRab stars 
to bend toward higher 04i values and to merge in the 
upper tail (see Clement & Shelton 1999); 

2. in NGC 6934, the RRab variables with 0.49 d< 
P <0.56 d seem to belong to the lower sequence and 
the others to the central one; 

3. in NGC 6362 there are some RRab stars with 
P > 0.60 d belonging to the upper tail and other RRab 
belonging to the central tail. 



Therefore, it looks evident that only the general trend of 
a large sample can be considered the fingerprint of the 
metallicity of the cluster, not the Fourier parameters ob- 
tained from a few stars. 

In the OGLE database the sample is three times more 
numerous and the separation between the tails looks bet- 
ter defined than in the Clement & Shelton (1999, their 
Fig. 8) sample. The RRab stars in oj Cen describe a sharp 
central tail (Poretti 2001, in preparation); a comparison 
with Fig. [ij supports the hypothesis that the three OGLE 
stars at P '^0.80 d belong to this tail. 

The presence of the tails in the OGLE database in- 
dicates an additional influence apart from metallicity on 
the light curve shape. We can conjecture that the shapes 
of the light curves are related to something more general, 
such as the evolutionary phase, or a combination of differ- 
ent physical conditions, since no obvious metallicity-based 
segregation is expected for field RR Lyr stars. 



5.3. RRc stars 

The slope in the (psi—P plot is steep. The relationships 
given by Simon & Clement (1993) allowed us to calculate 
some physical parameters. Morgan et al. (1998b) found 
a very homogeneous group, with a few stars showing un- 
usual metallicity. We subdivided the RRc sample into two 
subsets, the first (21 stars) with 0.25< P <0.30 d and the 
second (17 stars) with 0.30< P <0.35 d. There is no dif- 
ference in luminosity {logL— 1.65±0.01 and 1.67±0.01, re- 
spectively) and in the relative helium abundance (y=0.29 
everywhere) between the two subsets. On the other hand, 
the RRc stars with the shorter period seem more massive 
(mean mass 0.57±0.02Mq vs. 0.50±0.02 Mq) and hotter 
(mean temperature r=7460±10 K vs. r=7380±10 K), 
but the differences can be considered marginal. However, 
this tendency is confirmed when comparing the two RRc 
variables having P <0.25 d (0.56 and 0.57 Mq, 7726 K 
and 7649 K) with the two RRc stars having P >0.35 d 
(0.58 and 0.78 Mq, 7177 K and 7047 K). 

When considering different globular clusters we note 
that metallicity is driving the differences in the physi- 
cal parameters (the mean masses and mean luminosities 
increase and the mean temperatures and helium abun- 
dances decrease with decreasing cluster metal abundance, 
Clement & Shelton 1997). In the OGLE database and in 
the same globular cluster (i.e. under the same metallic- 
ity), the RRc variables with shorter periods are the more 
massive and hotter stars (see V2 and VI 1 in M55, Olech 
et al. 1999, or V14 and V28 in NGC6362, Olech et. al. 
2000) 

This tendency needs to be better fixed by a more ro- 
bust observational sample, which consider also RRc stars 
with P < 0.25 d (i.e. the extension of the RRc class we 
introduced here) and RRc stars with P > 0.35 d. 
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Table 2. Distribution of the cluster RRab light curves in the different segments (common part with P < 0.50 d and 
three tails) of the 04i— P diagrams. 
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-1.29 


30 


13 
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14 


M9 


-1.72 


6 




6 






M55 


-1.81 


6 


1 


3 
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5.4. The comparison with the other field RR Lyr stars 

Simon & Teays (1982) made the first Fourier decompo- 
sition of field RR Lyr stars, including both RRab- and 
RRc-types (70 stars). Subsequent papers did not further 
discuss their conclusions. As regards RRab stars, Simon & 
Teays claimed evidence for a trend, a sort of Cepheid-like 
progression, for P > 0.55 d. This progression is masked by 
a large scatter, not observed in the OGLE sample. As a 
consequence, these authors missed the opportunity to note 
that the progression originated earlier, around P ~0.45 d 
and to note the different tails. The reason for such a scat- 
ter is unclear, since observational errors or uncertainties 
on the parameters are not able to explain it fully. 

Simon & Teays also emphasized the unusual case con- 
stituted by XZ Get (P=0.8231 d), by far their longest- 
period star (gap from 0.75 to 0.82 d). The OGLE sample 
fills such a gap only partially (four stars clustering around 
0.77 d) and individuates a longer-period variable, MM7- 
B_V12 (P=0.8812 d). Gonsidering the whole sample of 
long-period RRab stars, the general pattern seems to be 
a continuous progression, not interrupted by any disconti- 
nuity, as suggested by Simon & Teays (1982). We also note 
that these long-period stars are not found in the globular 
clusters we considered above: in Figjl^ the longest period 
is 0.72 d. As a matter of fact, they seem to characterize 
metal-rich globular clusters, as NGG 6388 and NGG 6441 
(Pritzl et al. 2000). 

6. Conclusions 

The Fourier coefficients with their error bars are reported 
in a machine readable form in Tab. 3 (HADS), Tab. 4 
(RRc) and Tab. 5 (RRab). Moreover, an atlas of the light 
curves, least-squares fit and Fourier parameters is avail- 
able as a ps-file from the author. 

The step-by-step analysis we made allowed us to han- 
dle a sample characterized in an unambiguous way: no 
double-mode or Blazhko pulsator, no under- or over- 
evaluation of the significance of the Fourier components. 
Gonsequently, clear outputs have been obtained. 

Th e utility of the Fourier decomposition as a classifi- 



and amplitude parameters in well-defined loci for each 
class of variables. The use of this technique as a link be- 
tween observable quantities and stellar physical parame- 
ters is demonstrated by the enhanced separations within 
the same class, as examplified by the three tails in the 
RRab phase parameter plots. This way of performing the 
analysis when managing limited samples should give some 
warnings about the automated analysis of very large sam- 
ples. A good trade-off could be the introduction of sophis- 
ticated checks at intermediate steps. 

In a more general way, the exploitation of the large 
number of light curves supplied by the microlensing 
projects should be carried out in a careful way, since it 
can really provide us with the possibility to sound the 
stellar interiors. 
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